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Mixed Heisenberg Chains. I. The Ground State Problem. 
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We consider a mechanism for competing interactions in alternating Heisenberg spin chains due to 
the formation of local spin-singlet pairs. The competition of spin-1 and spin-0 states reveals hidden 
Ising symmetry of such alternating chains. 
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Introduction 

During the last few years mixed quantum spin chains 
have attracted some interest of theorists. Exactly solv- 
able versions with sophisticated Hamiltonians have been 
studied via Bethe ansatz Very recently numeri- 

cal methods |||| and matrix-product techniques have 
been applied to these spin systems. Different kinds of 
alternating chains with XXZ-like interactions have been 
investigated by using finite-size calculations and confor- 
mal invariance [Q. The subject of interest in this paper 
are alternating spin chains, in which each second site of 
the chain is considered as a compound, a kind of dumb- 
bell configuration. Two dumbbell spins 1/2 interact with 
each other via the Heisenberg interaction with a coupling 
constant Jo, either ferromagnetically, or antiferromagnet- 
ically. Each first site of the chain is occupied (A) by a 
usual quantum spin (we shall consider the cases of spin 
1/2, 1, 3/2, 2), or (B) by a compound spin too (see Fig. 
1). This spin is supposed to be antiferromagnetically cou- 
pled (coupling constant J\ < 0) to the spins of nearest 
dumbbells. In spite of short range Heisenberg interac- 
tions, variations of Jo may result in first order transitions 
at zero temperature. 

Realizations of such one-dimensional chains are shown 
in Fig. l(a and b). The Hamiltonian in case A can be 
written as 



and 



<p,r> 

-Jo J2 ^i) • o-(r 2 ), (1) 



<ri,r 2 > 

whereas in case B it becomes 



Hib) = E ( s (^) + s (^)) • Wn) + <r(ra)) 

<p,r> 

-Jo E <Pi) ■ <P2) ~ Jo E o-(n) • o-(r 2 ). (2) 



<Ti ,r 2 > 



<Pl,P2> 

Except for one special case which is described by model 
B, we concentrate our efforts on model A. The Hamilto- 
nians in (|[Q) can also be represented as follows: 

H (a.b) =Hi + Hi = _ Jx j2 s(p) ■ S(r), (3) 

<p,r> 



H, 



(b) 



^ = -^E s2 M' ( 4 ) 

r 

ijo'^S 2 ( P )-ijoE s2 W- ( 5 ) 



When making transformations from Hamiltonians (|l|j^) 
to we have skipped irrelevant constant terms. Co- 

ordinates of spins in dumbbells, r\ and r 2 , are replaced 
by a single co-ordinate r (in model B p\ and p 2 are also 
replaced by their common p). Note that model B trans- 
forms into model A when Jk — ► oo. 



a) 




FIG. 1 Shown are chain fragments of length k = 4. Model 
A: Spins s (= 1/2, 1, 3/2, 2) occupy sites p, spins a are ar- 
ranged in dumbbells. Model B: Spins s are in dumbbells, 
too, which are shown orthogonal to the tr-dumbbells. 

The two spins, <r(r{) and <x(r 2 ), are incorporated into 
the compound spin: S(r) = <T(rx) + o"(r 2 ), which is either 
0, or 1. This reveals a hidden Ising symmetry of the orig- 
inal Heisenberg models ([!]) and (g). In fact, Hamiltonian 
Hi does not generate any transitions between the total 
spin states and 1 of any compound spin. Thus we can 
introduce spin-0 states on (some) r-sites which are kind of 
intrinsic "defects" . Governed by the Jo-terms, these "de- 
fects" regulate a separation of the original chain into an 
ensemble of finite chain fragments decoupled from each 
other. Their structure can be defined as follows: A chain 
fragment of length k (k > 1) consists of k + 1 spins s and 
k spins 1. The spins of these two groups alternate with 
each other. A chain fragment can be formally described 
as (s, l) k s. Chain fragments are isolated from each other 
by spins 0. 
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For convenience, we enumerate the lattice sites r by 
integer numbers, then half-integers are reserved for sites 
P- 

All energies are suitably measured in units of J7i which 
is supposed to be negative. Thus, we set J\ = —1. Be- 
low we shall determine the phase diagram of the chains 
with s — 1/2, 1, 3/2, 2 depending on the parameter Jo 
in model A, and both, j7o and J"q, in model B. 

Our analysis includes elements of a rigorous analytical 
approach, a linear programming method, and numerical 
methods. Here we use the density matrix renormalization 
group (DMRG) method ||,|9| which is most appropriate 
to our problem. 

A detailed description of the DMRG algorithm which 
we have used to compute ground state energies of finite 
open chains can be found in Appendix A. It was imple- 
mented in C++ and ran on a SUN UltraSparc 2 worksta- 
tion with two 167 MHz processors and 256 MB memory. 
In order to achieve the desired accuracy we kept up to 
N = 100 block states during each DMRG step. The 
whole project consumed about 400 hours of CPU time. 
No parallelization was used. 



The ground state problem. Model A 

For convenience, we include into the definition of a 
chain fragment of length fc a spin-0 state, say, from its 
right. Then a chain fragment of length fc is represented 
by (s, l) fe (s, 0). This classification needs the "empty" 
chain fragments to be included: Any of these "empty" 
fragments is spin s with spin 0, attached from its right, 
i.e. (s,0). Conventionally, a nearest spin from the left of 
any chain fragment is also 0, but it is incorporated into 
the nearest-from-left chain fragment. 

Let us suppose that we have succeeded to determine 
the ground state energies of Hamiltonian Hi for all finite 
chain fragments, i.e. {eo, ei, £2, ■ ■ ■ , £fe> ■ ■ ■}• Then, the 
contribution of H^ will be — fcj7o for the chain fragment 
of length k. The ground state energy (per compound 
spin) of the chain, consisting of Nq "empty" chain frag- 
ments, N% chain fragments of length 1, . . . , Nk chain 
fragments of length k, etc. is 

£ g .s = y^(efc - kJo)wk- (6) 

fc>0 

In (H) we introduced "probabilities" Wk = Nk /A/", k > 
and J\f being the total number of compound spins. Ex- 
pressing Af in terms of the numbers Nk of the various 
chain fragments is equivalent to the constraint 



i = ^T(k + i)wk 



(7) 



fc>0 



imposed on the set of "probabilities" . Of course, all Wk 
are non-negative. Eqs. (^|)-(0) constitute a linear pro- 
gramming problem which demands that extrema of the 



energy (j^) should be searched at the vertices of the poly- 
gon defined by (^). This method has been efficiently ap- 
plied to the problem of competative interactions, leading 
to complex modulated structures 10-1^]. For this par- 
ticular problem, it can be easily proven that each vertex 
is characterized by only one non-zero "probability" value. 
For example, the vertex with wq = 1 corresponds to the 
perfect structure with the periodicity element (s,0), its 
energy being £ — (we denote this regular spin configu- 
ration by (0)); w\ — 1/2 (a configuration conventionally 
denoted by (1)) corresponds to the periodicity element 
(s, 1, s, 0) with energy £ = (ei — Jo)/2. For the (k) state, 
Nk = N '/ (fc + 1) and the periodicity element can be rep- 
resented as ((s, l) fc , s, 0). The energy of this spin config- 
uration is given by 



£k = (e fc - kj )/(k + 1). 



(8) 



Numerical methods which are outlined in Appendix A 
allow us to analyze characteristic first order transitions. 
They happen at zero temperature and are controlled by 
Jo- For s = 1/2 the set of energies {e^} is given in Table 
I (left column). According to the Lieb-Mattis theorem, 
a ferrimagnetic ground state would be realized with a 
total spin fc/2, if we dealt with a periodic alternating 
chain, consisting of k spins 1/2 and k spins 1, coupled 
antifcrromagnctically. In our case, a chain fragment with 
k spins 1 and k + 1 spins 1/2 will exhibit the total spin 
Sk = (fc — l)/2 (fc > 1) in the ground state. Low lying 
excitations are asymmetric. Their hierarchy is as follows: 
the lowest excitations are triplets which correspond to 
AS = -1, or the total spin (fc - 3)/2 (fc > 3). The 
next lowest excitations are triplets too, but with AS = 
+1 (the total spin (fc + l)/2). The singlet excitations 
are lying above both triplet excitations. Only the triplet 
excitations with AS = — 1 will give rise to a gapless mode 
in the limit fc — > oo. 



e. (m) 



FIG. 2 Shown are a few first points of ei n t(w) (model A, 
s = 1/2). 

Comparison of the ground state energy expressions (|^) 
at various vertices clearly shows that three configurations 
are competitive in the global ground state: This is (0) for 
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Jo <— 2, which changes to (1) for —2 < Jo <2e QO + 2 w 
-0.910. Beyond this region, i.e. for -0.910 < Jo, the 
(oo) state becomes energetically favourable. The numer- 
ical data for the set {ek} can be represented as 



— ke c 



e 



e int (fc), 



(9) 



where eoo ~ —1.45412 is the energy per element (s, 1, s) 
of the perfect periodic spin structure (oo) = (s, 1)°°, s = 
1/2. The energy due to the open ends is e ~ —0.45352, 
and the remaining part, ej n t(fc), can be interpreted as the 
interaction between the chain fragment ends which ap- 
parently turns to zero for k — > oo. This function e- m t(k) 
plays an important role in establishing the succession of 
phase transitions. In Appendix B we perform a rigor- 
ous analysis, according to which a succession is given 
by a broken line, which is concave upwards and envelops 
eint(fc) from below. This broken line includes the points 
corresponding to periodic structures with shortest chain 
fragment (fc = in model A) and infinitely long chain 
fragments, k — + oo. 





s — 


1/2 


s — 


3/2 


s — 
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ft 








ei„t(fe) 




ei„t(fe) 








0.453733 





0.150038 





0.111298 


1 


-2 


-0.092173 


-4 


0.011951 


-5 


0.005167 


2 


-3.3815016 


-0.019581 


-7.8727492 


0.001114 


-9.8987531 


0.000282 


3 


-4.8191715 


-0.003157 


-11.735772 


0.000003 


-14.792914 


-0.000009 


4 


-6.2705355 


-0.000426 


-15.597742 


-0.000054 


-19.686790 


-0.000017 


5 


-7.7242539 


-0.000051 


-19.459627 


-0.000026 


-24.580654 


-0.000012 


6 


-9.1783024 


-0.000005 


-23.321522 


-0.000009 


-29.474518 


-0.000007 


7 


-10.632391 


0.000000 


-27.183429 


-0.000003 


-34.368384 


-0.000004 


8 


-12.086485 


0.000000 


-31.045338 


0.000000 


-39.262251 


-0.000003 


9 


-13.540579 


0.000000 


-34.907251 


0.000000 


-44.156119 


-0.000002 


10 


-14.994673 


0.000000 


-38.769163 


0.000000 


-49.049987 


-0.000001 



TABLE 1 The ground state energies of H\ and the function 
eint(fc) for the chain fragments, consisting of k spins 1 and 
k + 1 spins s, 1/2, 3/2 and 2 in three successive columns. 



ft 


e2fc 


e tat (2fc) 


€2fc + l 


e int (2fc+l) 








or € fc 


or ei n t(fe) 








-1.208 





0.193484 


1 


-2 


-0.405032 


-3 


-0.003548 


2 


-4.6457513 


-0.247815 


-5.8302125 


-0.030792 


3 


-7.3702750 


-0.169371 


-8.6345320 


-0.032144 


4 


-10.124637 


-0.120765 


-11.432932 


-0.027575 


5 


-12.894560 


-0.087720 


-14.230359 


-0.022035 


6 


-15.674010 


-0.064202 


-17.028266 


-0.016973 


7 


-18.459853 


-0.047076 


-19.827036 


-0.012775 


8 


-21.250218 


-0.034473 


-22.626683 


-0.009455 


9 


-24.043879 


-0.025167 


-25.427099 


-0.006902 


10 


-26.839978 


-0.018298 


-28.228140 


-0.004975 


11 


-29.637889 


-0.013240 


-31.029675 


-0.003542 


12 


-32.437147 


-0.009530 


-33.831576 


-0.002474 


13 


-35.237402 


-0.006818 


-36.633771 


-0.001701 


14 


-38.038394 


-0.004842 


-39.436175 


-0.001138 


15 


-40.839927 


-0.003407 


-42.238735 


-0.000730 



TABLE 2 The ground state energies of Hi and the function 
eint(fc) for the chain fragments, consisting of spins 1. For 
model B, the data of computations are given in the left (right) 
column for chain fragments of even (odd) total lengths, 2k 
(2k + 1). For model A, only the numerical data of the right 
column should be taken. 



For s = 1/2, ei n t(fc) is shown in Fig. 2. Thus, the en- 
ergetically favourable configurations can be (0), (1) and 
(oo). 

For a = l, the function ei nt (fc) is the upper curve in 
Fig. 3 related to odd integers m = 2k + 1. According to 
Appendix B we have successive phase transitions 

(0) - <1> - (2) - <3> - <oo) 

at J (0,1) = -3, J (1 ' 2) = -2.660425, J (2 ' 3) = -2.582746, 
J (3 ' oo) = -2.577340, which are determined from Eq. (§) 
and data of Table 2 (right column) . In accordance with 
the Lieb-Mattis theorem, the value of the spin of the 
ground state of a chain fragment is Sk = 1, independent 
of length. 

For s = 3/2, the total spin in the ground state is Sk = 
(3 + k)/2, whereas for s = 2, Sk — k + 2. A succession of 
phase transitions can be easily identified by making use 
of Table 1 (middle and right columns). It is the same in 
these two cases: 

(Q) - <1) - <2) - (3) - (4) - <oo). 

The last two transitions occur at Jq 3 ' 4 ^ and Jq 4 ' 00 -* whose 

(2 3) 

values slightly differ from J V ' '. 

Note that a transformation (0) — ► (oo) via a few inter- 
mediate first order transitions proceeds due to quantum 
effects. In fact, if Hi is confined to a pure Ising form, 
then ei n t(fc) = 0. The system undergoes a direct transi- 
tion (0) — > (oo) at Jo = — 2s. A transition point is highly 
degenerate: Chain fragments of any length are allowed as 
well as any sequence of them. 

In contrary to this, a spin- wave approach "overesti- 
mates" quantum effects for s > 3/2: e- m t(k) of the spin- 
wave approach may not differ much from the numerical 
data, but for s > 3/2 it always decays monotonically with 
k and thus would lead always to an infinite set of first 
order transitions (0) — ► (1) (k) — ► ... — * (oo). 

Exceptional is the case s = 1/2 for which the spin- 
wave approach exhibits a global minimum of ei n t(fc) at 
k = 1. This results in the same succession of transitions 
(0) — * (1) — ► (oo), as obtained within the exact numerical 
scheme. 



The ground state problem. Model B 

The subject of this Section is model B, where a and s 
are both spins 1/2, thus the equivalent model described 
by Hamiltonian <j3j> deals with chain fragments consist- 
ing of compound spins 1 only. The ground state energies 
and excitations of finite spin-1 chains described by Hamil- 
tonian Hi have been studied by T.Kennedy |l4| in the 
framework of the Lanczos method. 

It is more suitable to enumerate the chain by integers, 
say odd and even for p and r sites, respectively. Our 
consideration can be restricted to J ' > Jo ■ Relation (||) , 
which can be used for spin-1 chain fragments as well, is 
apparently not identical for chain fragments consisting of 
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even or odd number of sites. This does not concern e^, 
whose value is common for both types of chain fragments, 
but concerns cq JlH] and e; n t(A;). The total spin in the 
ground state is zero or one, it depends on the length of 
chain fragments, even or odd. 

The linear programming method will be used again 
to select those chain fragments which may be candi- 
dates to form the ground state. The following four 
types of fragments should be taken into consideration: 
(l 2fe+1 ,0) r , (l 2fc+1 ,0) p , (l 2fe ,0) r , (l 2fe ,0) p . We de- 
note their total numbers by N$ +1 ,N$ +1 ,N$,n!&\ 
respectively, and the corresponding "probabilities" by 



(r) 

°2k+l> 



(p) (r) (p) 
°2h+l> w 2k > W 2k . e 'g- 



lj 2k+l 



N^'/N. Here 



Af is the total number of sites of both types, r and p. 
As in the prior consideration, a spin-0 state is attached 
from the right to any finite spin-1 chain fragment. In- 
dices r and p specify the type of the rightmost site, which 
is occupied by spin-0. Certainly, this definition incorpo- 
rates "empty" chain fragments r and P into the scheme. 
Not all the numbers N mentioned above are independent, 
e.g., the leftmost sites of (l 2fc , 0) r and (l 2fe+1 , 0) p are both 
of the r-type. They evidently follow all those chain frag- 
ments, which have a rightmost site of type p. Thus, the 
following "conservation law" holds: 



(r) 



N, 



(p) 

2k+ 



1 ) = E^ ) = 



i 



-2fc 



2fc + 1 



(e 2 fe - kj ' -kJo). 



(13) 



For the second kind we get w^fe+i = l/(2fc + 2) and the 
energy reads 



1 



2fc+l 



2k 



:{£2k+i-{k + l)J^-kJ ). 



(14) 



We shall denote corresponding regular structures as 
((2k)) and ((2k + 1)) in accordance with the numbers 
of spin-1 sites in chain fragments. Note that in the sec- 
ond case the periodicity is 2k + 2, whereas in the first 
case it is 4fc + 2. Expressions (|l^-(|l4|) can be rewritten 
as 



f 



1 



771+1 



Jo + Jo 



Jo — Jo 



, (15) 



where 8 m — or 1 for m even or odd. Except for the last 
term in the r.h.s., Eq.(|l5|) has a form similar to Eq.(H). 

e ; _ fm) 



which yields 



2k ■ 



(10) 



2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 
• ••• •• m 



The total number of lattice sites M expressed in terms of 
and gives rise to the equation: 



1 = ^(2(2^ + 1) 



+ (2fc + 2)(m 



$H+«4&-i)> (n) 



k>0 



and the energy per site is 



_>/,- 



+ (£2k+i - (k + 1)X ~ kJo)w ( £ +1 



+ {e 2k +i - kj ' - (k + l)Jo)w. 



2k+l 

(P) 

2k+l ■ 



(12) 



In this equation we should use energies e m given in Table 
2. Note, that the set {£2fc+i} coincides with {e^} used in 
model A. 

As we have assumed J ' > Jo , the contribution of chain 
fragments (l 2fc+1 ,0) p to the energy is not competative 
with the one for the (l 2fc+1 ,0) r fragments, see Eq. (12). 
Thus, in the problem of finding the ground state energy, 
we must check the contribution of two sorts of vertices: 
The first kind is defined as = w$ = l/(4fc + 2) with 
energy 



FIG. 3 Shown are two subsets of ei n t(m) for m odd and 
even. The points for m = 2 and m — are situated far below 
at -0.4056 and -1.2086, respectively. ei nt (2fc + 1) of model B 
coincides with emt(fc) of model A. 

Shown in Fig. 3 is e- m t at J Q ' = Jo- If, however, J Q ' > 
Jo, then the subset {ei n t(w, even)} increases as compared 
with {e- m t(m, odd)}. Having Fig. 3 as a prerequisite, we 
can describe all possible transformations as a function of 
(Jo -Jo): 

• As far as ej nt (0) remains below e; nt (7), which is 
the absolute minimum of the {e m t(m, odd)} sub- 
set, our system only undergoes the ((0)) — > ((oo)) 
transition; 
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• If ei n t(0) raises above ei n t(7), but remains below 
the continuation of the straight line connecting the 
points at m — 5 and at m = 7, then two transitions 
((0» -» ((7)) -» ((oo» take place; 

• Next, if ei n t(0) raises above the continuation of the 
line connecting the points at m = 5 and m = 7, 
but remains below the continuation of the line con- 
necting the points at m = 3 and m = 5, then the 
three transitions ((0)) -> ((5)) -> ((7)) -> ((oo» 
take place. 

• The scheme can be continued. 

• Finally, at sufhciently large Jq — Jo, we get a 
maximal possible number of first order transitions: 
((0)> ((1)) - ((3)) - ((5)) -+ ((7)) - «oo». 

The phase diagram is shown in Fig. 4. 



«1» 




-3.1 -3 -2.9 -2.8 -2.7 -2.6 -2.5 

FIG. 4 The phase diagram for model B is shown in the 
{Jo, Jo) plane. A narrow area between ((5)} and ((oo)) be- 
longs to the ((7)} phase. The line, which separates the ((0))- 
phase from the others, changes its slope from —1 through 
-3/4, -2/3, -1/2 to 0. 



Discussion and Conclusions 

The model considered in this paper is probably a sim- 
plest representative of the family of Hcisenberg models 
which possess an intrinsic property of a hidden Ising sym- 
metry. Compound spins of our model form a regular 
sublattice within a one-dimensional chain: We have con- 
sidered only the case of alternating chains, however, any 
periodicity within a system of compound spins is allowed. 

We have restricted ourselves to the compound spin val- 
ues and 1. In principle, a more complex construction for 



compound spins can be used, i.e. instead of dumbbells 
three spins 1/2 may form a triangle. Then a compound 
spin is allowed to be 3/2 and 1/2. A treatment of this 
case is more difficult, because all possible modulations in 
a distribution of spins 3/2 and 1/2 should be considered 
within the ensemble of infinite chains, there are no spin-0 
breakers like in the models of this work. Another pos- 
sibility would be associated with a constituents of com- 
pound spins of a higher value than 1/2, say 1. Then, 
increasing Jo from large negative values to moderate neg- 
ative ones, we could pass through a few regimes, starting 
from a periodic structure of elements (s, 0) through a few 
structures whose periodicity elements are ((s, l) fe , s, 0) to 
(s,l)°°, then, most likely, a few intermediate structures 
(s,2) fe ,s, 1 will finally lead to a perfect (s,2)°° chain. 
However, the analysis of the phase diagram between per- 
fect (s, 1)°° and (s,2)°° structures is somewhat difficult 
because of the absence of spin-0 breakers. 

For dealing with these more complicated systems, one 
could use a spin-wave approach, although tedious for fi- 
nite systems, but well-defined. We have employed this 
approach in order to compare our "exact" numerical 
results with this approximate analytical scheme. For 
s = 1/2 these two schemes qualitatively lead to the same 
global minimum of e; nt (fc) at k = 1, thus exhibiting the 
identical succession of transitions. We have observed that 
for s > 3/2 eint(fe) is always a monotonically decreasing 
function. This means that the occurence of a negative 
minimum in ej n t(fc) should be due to non- linear terms 
of a spin-wave expansion. It is clearly seen from Table 
1 that this minimum is extremely small (~ 10~ 4 ei n t(0)) 
even at s = 3/2 or 2, but it exists and has an impor- 
tant influence on the succession of phase transitions vs 
Jo- As this minimum occurs at k = 4, which is a rather 
short length, and thus is well-controlled numerically, we 
are absolutely sure of the existence of the minima. The 
data obtained from the spin-wave approach are given in 
Table 3. Evidently they do not differ much from our 
"exact" numerical results, decreasing very rapidly with 
k. Exceptional is the s = 1 case, where ei nt (fc) decreases 





s — 


1/2 




/2 




k 


(■k 


eint(fe) 




eint(fe) 




eint(fe) 








0.370671 





0.184008 





0.131569 


1 


-2 


-0.192873 


-4 


0.012081 


-5 


0.004480 


2 


-3.2928932 


-0.049311 


-7.8377223 


0.002433 


-9.8768944 


0.000497 


3 


-4.6937567 


-0.013719 


-11.667585 


0.000644 


-14.750229 


0.000073 


4 


-6.1203597 


-0.003866 


-15.496108 


0.000194 


-19.623201 


0.000012 


5 


-7.5540280 


-0.001079 


-19.324313 


0.000063 


-24.496123 


0.000002 


6 


-8.9897014 


-0.000297 


-23.152429 


0.000022 


-29.369036 


0.000001 


7 


-10.425941 


-0.000081 


-26.980516 


0.000008 


-34.241947 


0.000000 


8 


-11.862337 


-0.000022 


-30.808595 


0.000003 


-39.114859 


0.000000 


9 


-13.298777 


-0.000006 


-34.636671 


0.000001 


-43.987770 


0.000000 


10 


-14.735228 


-0.000002 


-38.464745 


0.000000 


-48.860681 


0.000000 



TABLE 3 Spin-wave approximation: The ground state ener- 
gies of Hamiltonian H\ and the function ei nt (fc) for the chain 
fragments, consisting of k spins 1 and k + 1 spins s, 1/2, 3/2, 
and 2, are shown in three successive columns. 
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slowly, as within the "exact" numerical scheme, although 
monotonically. 

We suppose that this important property of ei n t(fe) 
to exhibit a minimum, more and more shallow, will be 
valid for larger s-values, too, but we cannot say definitely 
whether the minimum remains at finite fc-values or shifts 
to k — > oo. 

The next article in this series which will be published 
elsewhere is devoted to the thermodynamics of mixed 
Heisenberg chains. Certainly, all the transitions will be 
smeared out due to thermal fluctuations. However, the 
system must show big changes in physical properties, 
such as the specific heat and the magnetic susceptibility, 
if J7b is in a transitional area. In fact, as the chain frag- 
ments of finite lengths possess total nonzero spins, which 
do not interact with spins of the nearest chain fragments, 
they form a system of paramagnetic spins. Thus, the 
susceptibility will exhibit a Curie-like behaviour at low 
temperatures. A prefactor will be temperature depen- 
dent, too, because the values of paramagnetic spins and 
their concentrations strongly correlate with the length 
distribution of the chain fragments. 

We illustrate this by taking an example from the forth- 
coming article. Let us consider model A with s = 1/2 
and j7o ~ —2, i.e. where the chain fragments (s,l,s,0) 
and (s, 0) are only competitive. The former is practically 
in a singlet state, the latter represents a purely param- 
agnetic spin 1/2, their concentration in a lattice varies 
with T as 



1 1 + VI + 4w 



2 1 + Aw + yT + 4w 



= exp(Jo + 2)/T. 
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Each box denotes a single-spin Hilbert space with basis 
states numbered by the corresponding index i v in the 
lower row. For the moment ii, 13, and 15 are identical, as 
well as 12 and i^. The Hamiltonian for the chain (A.l) 
can be divided into the following contributions: 



H 
H 



BS2 I TJBS'2 



H S2 - +H S2 - +H Sl - 

12,32 ' 'I4J4 ' 13,33 

I TTSIS2 I TJ-S1S2 

1312,3332 i3»4,j334' 



(A.2) 



The first five terms are on-site contributions, the others 
couple neighbouring sites. At this stage, H B , H BS2 are 
identical to H Sl , H SlS2 , respectively. 

The first step of the algorithm now is to compute the 
lowest eigenvalue of the Hamiltonian (A2) and the cor- 
responding eigenstate by using the Lanczos method. The 
eigenvalue is directly the ground state energy 62 of the 
fragment (Si, S2) 2 Si and the corresponding eigenstate 
4>iii2i3iiia serves as the 'target state' of the subsequent 
DMRG step. From this target state we calculate the 
density matrix for the combined i\i2 Hilbert space 



Phi2,jih 



E • 

i3,ii,i5 



(A.3) 



Clearly, the eigenvectors of p to the largest eigenval- 
ues give the most important contributions to the target 
state <bi,i^i,i.i,. The idea of the DMRG method is to 



reduce the dimension of the i\i2 Hilbert space by pro- 
jecting all operators onto the eigenstates of p belonging 
to the N largest eigenvalues. Let U be the rectangular 
matrix that contains these TV normalized eigenvectors as 
columns |l6| . This matrix U allows us to combine the i\ 
and 12 Hilbert spaces into a new single Hilbert space i\ 
while limiting its dimension to at most N. All operators 
A acting on ziZ2 are transformed into operators A' acting 



on i\ via 



A' = UAU^ 



We then arrive at the situation 



Si S2 Si S2 Si 



13 



(A.4) 



(A.5) 



Appendix A. 

To calculate the ground state energies of finite chain 
fragments (Si, S2) k Si we used a slightly modified version 
of the well-known density matrix renormalization group 
(DMRG) method introduced by White §,|. Here we 
give a brief sketch of this method. 

Let us consider the following finite fragment, which 
consists of 3 spin-Si sites and 2 spin-S2 sites: 



Si S2 Si S2 Si 

h «2 h u h 



(A.l) 



Note that i' 5 is the reflected version of i[. The length of 
the chain can now be increased by inserting an (Si, S2) 
pair, yielding 



Si S2 Si 



Si S2 Si 



(A.6) 



where we have renamed to i' 2 . The Hamiltonian for 
the chain (A.6) again has the structure (|A.2j), but with 



interchanged roles of Si and S2 and new operators 



Hi 



U [H B (g)l + l<E>H S2 + H BS2 } W 



(A.7) 



G 



provided the original interaction is given by 



Taking into account inequalities (B.4), one simply ob 



H 



SiS 2 



(A. 



This completes the DMRG step. The Lanczos met hod is 
now applied to the new full Hamiltonian for chain ( |A.6| ), 
the lowest eigenvalue gives the energy 63 of the chain 
fragment (Si, S^) 3 ^ and the corresponding eigenvector 
serves as the target state for the subsequent DMRG step. 

The algorithm described above is a straightforward 
generalization of the standard 'infinite system method' 
of White H to alternating symmetric spin chains. It 
differs from the original algorithm in two points: 1. It 
contains an additional spin 13 in the center of the chain, 
and 2. the role of Si and S2 has to be interchanged 
at every DMRG step. We did not use the extended 'fi- 
nite system method', as agreement with competing pure 
Lanczos calculations was already satisfactory and could 
be systematically increased by using larger N. 



Appendix B. 

In this Appendix we show how a succession of first 
order transitions takes place. We start by noting that 
the (0) structure is optimal when j7o has a large negative 
value and the (oo)-phase is realized at large positive Jo- 

Let us define the J7o parameter as Jg m ' n \ when two 
phases, (m) and (n), are at equilibrium, i.e., £ m = £ n . 
Then we can employ representation (^) to determine the 
following equations: 



= (k + l)e lnt (fc - 1) - ke int (k), 

= (k + 2)e int (fc) - (k + l)e int (fc + 1). 
From Eqs.(pT|)-(^2|) we obtain 



(B.l) 
(B.2) 



„{k,k+l) „y 

(k + l)(e int (fc + 1) - 2e int (fc) + e int (fc - 1)). 



r (k-l,k) 



Thus, we come to the first conclusion: If ei n t(fc) is con- 
cave upwards at any k, we have a full succession of tran- 
sitions: (0) -> (1) -»• (2) -»• . . . -> (00). 

Suppose now that the condition for e; n t(fe) does not 
hold. For instance, if for < j < m e-^ik + j) satisfies 
inequalities: 

j 

eint(fc+i) > e int {k) -\ {e int (k + m) - e lnt (fc)), (B.4) 

m 



then a generalization of Eqs. (B.1)-(B.2) is 

eint(fc) 



r (k,k+j) _ 



eo 



k + 1 



tains j( k ' k+ ^ > j{ k ^ k + m ) _ xhese last inequalities show 
that all intermediate phases whose "energies" e- m t(£) are 
above the enveloping line, cannot be realized as ground 
states at any value of J . An important consequence is 
that if ei nt achieves an absolute minimum at some fc, this 
results in a direct transition (k) — * (00). A convenient 
expression for the critical value J7"g fe '°°'' can be derived 
from (B.5): 



^(fc.00) = j(k-l,k) + {k + 1)(e . nt(fc _ 1} _ e . nt(fe)) . (B _ 6) 
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